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Abstract The effective quark interaction in a relativistically covariant constituent quark model based on the 
Salpeter equation is supplemented by an extra phenomenological flavour dependent force in order to account 
for some discrepancies mainly in the description of excited negative parity A resonances. Simultaneously 
an improved description of some other features of the light-flavoured baryon mass spectrum and of some 
electromagnetic form factors is obtained. 



PACS. 11. 10. St Bound and unstable states; Bethe-Salpeter equations 
- 13.40. Gp Electromagnetic form factors 



12.39.Ki Relativistic quark model 



1 Introduction 



The description of the hadronic excitation spectrum re- 
mains a major challenge in strong interaction theory. In 
spite of recent progress in unquenched lattice QCD ac- 
cess to excited states is still very limited [UH]. Therefore 
it seems worthwhile to improve upon constituent quark 
model descriptions, which in view of the light quark 
masses (even taken as effective constituent masses) have to 
be formulated in terms of relativistically covariant equa- 
tions of motion. About a decade ago we formulated such 
a quark model for baryons, see [SHS] on the basis of an in- 
stantaneous formulation of the Bethe-Salpeter equation. 
In this model the quark interactions reflect a string-like 
description of quark confinement through a confinement 
potential rising linearly with interquark distances as well 
as a spin-flavour dependent interaction on the basis of in- 
stanton effects, which explains the major spin-dependent 
splittings in the baryon spectrum. 

Such a model description should offer an efficient de- 
scription of masses (resonance positions) , static properties 
such as magnetic moments, charge radii, electroweak am- 
plitudes (form factors and helicity amplitudes) with only 
a few model parameters. As such they also offer a frame- 
work which can be used to judge in how far certain fea- 
tures could be considered to be exotic. This concerns e.g. 
phenomenological evidence for states with properties that 
can not be accounted for in terms of excitations of quark 
degrees of freedom, as is at the heart of any constituent 
quark model, but instead requires additional degrees of 
freedom as e.g. reflected by hadronic interactions. 
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A satisfactory description of the major features in the 
light-fiavoured baryonic mass spectrum could indeed be 
obtained. These include 

— the linear Regge trajectories with an universal slope 
for all flavours including states up to total angular mo- 
menta of J = ^ and excitation energies up to 3 GeV, 
see nig ; 

— the position of the Roper-resonance and three other 
positive parity excited nucleon states well below all 
other states of this kind. These can be largely ac- 
counted for by the instanton-induced force, the 
strength of which was chosen to reproduce the ground 
state N - A splitting [1[S]; 

— a plethora of electroweak properties which can be ex- 
plained without introducing any additional parame- 
ters, see [5HH]- 

Nevertheless some specific discrepancies remain; most 
prominent are: 

— the conspicuously low position as well as the decay 
properties of the negative parity Ai - (1405) resonance; 

2 

The calculated mass of this state exceeds the experi- 
mental value by more than 100 MeV; 

— there is experimental evidence [HI for excited negative 
parity Z\-resonances well below 2 GeV which can not 
be accounted for by the quark model mentioned above, 
see fig.m nor by any other constituent quark model we 
are aware of; 

— The mass of the positive parity Z\3 + (1600) resonance, 

2 

see also fig. [H the low value of which with respect to 
other excited states of this kind can not be traced back 
to instanton-induced effects, since these are absent for 
flavour symmetric states. 
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Fig. 1. Discrepancies in the A mass spectrum: The left part of 
each column represents the results obtained in model ^ of [3] 
in comparison with experimental data from the Particle Data 
Group [9] (right side of each column) , where lines are the res- 
onance position (mass) with the mass uncertainty represented 
by a shaded box and the rating of [9] indicated by stars. J 
and TT denote total angular momentum and parity, respectively. 
Small differences with respect to the results from fig. 3 of [4] 
are due to the fact that we obtained increased numerical ac- 
curacy by diagonalising the resulting Salpeter Hamiltonian in 
larger model spaces, see section for details. 



We therefore want to explore whether these deficiencies 
are inherent to the constituent quark model itself or can 
be overcome by the introduction of an additional quark in- 
teraction which improves upon the issues mentioned above 
without deteriorating the excellent description of the ma- 
jority of the other states. In view of the fact that the dis- 
crepancies mainly affect the Z\-spectrum, this additional 
interaction is likely to be flavour dependent. An obvious 
candidate in this respect would be a single pseudoscalar 
meson exchange potential as has been used as a basis of 
an effective spin-flavour dependent quark interaction very 
successfully by the Graz-group [TOHlT] . 



In the present paper we shall investigate various imple- 
mentations of the coordinate (or momentum) dependence 
of such interactions. The paper is organised as follows: 
After a brief recapitulation of the ingredients and basic 
equations of our Bethe-Salpeter model (for more details 
see [3]) in section [5J we discuss in section |3] the form and 
the parameters of the effective quark interactions used in 
this paper. Section 2] contains the results and a discussion 
of the baryon mass spectra in comparison to the results 
obtained before [HI^. In section [5] we present some re- 
sults on ground state form factors before concluding in 
section El 



2 Bethe-Salpeter model 

2.1 Bound state Bethe-Salpeter amplitudes 

The basic quantity describing three-quark bound states is 
the Bethe-Salpeter amplitude x defined in position space 
through 



- {0\T^a^{xl)^a,{x2)^aAx3)\P) , 



(1) 



where T is the time ordering operator, |P) represents the 
bound-state with total 4 momentum = of a baryon 
with mass M , |0) is the physical vacuum and ^P'a- (x^) de- 
notes single quark-field operators with multi-indices in 
Dirac, colour and flavour space. Because of translational 
invariance below we shall exclusively use relative Jacobi 
coordinates , in momentum space. The Fourier trans- 
form of the Bethe-Salpeter amplitudes: Xp 010203 (-P^iPr?) 
are then determined by the homogeneous 
Bethe-Salpeter equation compactly written as 



Xp = -iGo,P [K^^^ + K'^> ] XP , (2) 

where Kp represents the irreducible 3-quark-kernel and 
where Kp is defined by 



'(2) 



k=l 



xK 



(2) 



iPi.,pl)®isiy\i-Pv.) (3) 



in terms of the irreducible two-body interaction kernel 

(2) 

> each quark pair labeled by the odd-particle 

index k. Furthermore Go.p is the free 3-quark fermion 
propagator defined as 

xSUi+Pi + ^)®SUi-Pi + ^y^SUi-Pv) (4) 
in terms of full single quark propagators Sp . 



2.2 Model assumptions 

In view of the fact that the interaction kernels and the 
propagators are sums of infinitely many Feynman dia- 
grams, in order to arrive at a tractable model we make 
the following assumptions, mainly with the goal to stay 
in close contact with the quite successful non-relativistic 
constituent quark model: 

— The full propagators Sp are replaced by Feynman pro- 
pagators of the free form 



Ship) 



le 



(5) 
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tacitly assuming that at least some part of the self- 
energy can effectively be subsumed in an effective con- 
stituent quark mass mi which then is a parameter of 
the model. 

Obviously this does not account for confinement: This 
is assumed to be implemented in the form of an in- 
stantaneous interaction kernel which in the rest frame 
of the baryon is described by an unretarded potential 

y(3) . 

= ^^'Up«,p,,;p?',p^') ■ (6) 

Likewise we assume that two-quark interaction kernels 
in the rest frame of the baryon are described by 2-body 
potentials V^'^^ : 



P={MM) 



(7) 



Then, with a perturbative elimination of retardation ef- 
fects, which arise due to the genuine two-body interac- 
tions, sec [3 for details, one can derive an equation for 
the (projected) Salpeter amplitude 



. f dp^l dp° 



where 



with 



Here 



± A-{pi) ® A-{p2) A-{p3) 



(p) ± H,n{p) 



(9) 
(10) 

(11) 



are projection operators on positive and negative energy 
states and Vf projects on quark flavour / . Furthermore 
the quark energy is given by uj(p) = ^|pp + and 

HM:^l°il-p + m) (12) 

is the Dirac Hamilton operator. As shown in detail in [3] 
this equation can be written in the form of an eigenvalue 
problem 



erator is given by 

[nKi] {Pi,Pv) = noiPi,Pr,)KiPi.Pv) 

d'pi' d'p,' 



-l-yl+(P5,p^)7°«)7°®7"' 



(27r)3 (27r)3 

^^'^(p?,Pr,;P«',P^')<(p«',P^') 

+ ^_(P5,P^)7"®7%1 J 

V^^\p^;p^')®t1>ii{Pi.Pv) (14) 
-|- corresponding quark interations (23) and (31) , 

where "Hq denotes the free three-quark Hamilton opera- 
tor as a sum of the corresponding single particle Dirac 
Hamilton operators. 

The eigenvalue problem of Ec^. (fT3|) is solved by nu- 
merical diagonalisation in a large but finite basis of oscil- 
lator states up to an oscillator quantum number iV,nax ■ In 
previous calculations [4l[5] at least iVmax = 12 was used. 
All the results in the present paper were obtained with 
at least iVmax — 18 , which, although computer time con- 
suming, has the advantage that for all states the indepen- 
dence of the numerical results on the oscillator functions 
length scale in some scaling window could be warranted 
and that all could be calculated with a universal value 
for this length scale. This is a technical advantage when 
calculating electroweak amplitudes. 



3 Model Interactions 

Below we specify the interaction potentials V^^^) g^-^d V'^'^^ 
used in Eq. (ITil) . These include a confinement potential, 
the instanton induced two-quark interaction as has been 
used before [HIS] and the new phenomenological potential 
inspired by pseudoscalar meson exchange. 



3.1 Confinement 

Confinement is implemented by subjecting the quarks to 
a potential which rises linearly with interquark distances, 
supplemented by an appropriate three particle Dirac 
structure F. The potential contains two parameters: the 
off-set a and the slope h and is assumed to be of the fol- 
lowing form in coordinate space 



i<3 



(15) 



where Fo and Fg are suitably chosen Dirac structures. 
Alternatively we can consider the linear potential to be 
treated as a two-body kernel as will be used below. 



M 



(13) 3.2 Instanton induced interaction 



for the projected Salpeter amplitude, where the eigenval- 
ues are the baryon masses M . The Salpeter Hamilton op- 



Instanton effects leads to an effective quark-quark inter- 
action, which for quark pairs in baryons can be written in 
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coordinate space as 

Vlf^ixi,x2;x[x2) = Vl^{xi - X2) 



and 



III v-^ii -^2, ■^i-i^2) — "m 

X 5{x1 - x9,) S^*\xi ~ x[) S^*\x2 ~ x'2) (16) 



with 



®{9nn V^{nn) + gns V^{ns)) , (17) 

where V^^^^q is a projector on spin-singlet states and 
'^^(71/2) projects on flavour-antisymmetric quark pairs 
with flavours /i and f2- Ahhough the two couphngs g„„ 
and gns are in principle determined by integrals over in- 
stanton densities, these are treated as free parameters 
here. As it stands this is a contact interaction, which 
for our purpose is regularised by replacing the coordinate 
space dependence by a Gaussian 



vxix) 



exp 



(18) 



The effective range parameter A is assumed to be flavour 
independent and enters as an additional parameter. 



3.3 An additional flavour dependent interaction 

The coupling of spin-i fermions to a flavour nonet of pseu- 
doscalar meson fields is given by an interaction Lagrange 
density 



(19) 



in the case of so-called pseudoscalar coupling and by 



a=0 



(20) 



for pseudovector coupling. Here "0 represents the quark 
fields with mass m and 0° the pseudoscalar meson fields 
with mass fia where the fiavour index 
ivr°, K'^ . The fiavour dependence is represented by the 
usual Gell-Mann matrices A°, a — 1, . . . , 8 ; A° is propor- 
tional to the identity operator in flavour space normalised 
to Tr((A°)2) = 2. 

A standard application of the Feynman rules, see e.g. 
[15] then leads to the second order scattering-matrix ele- 
ment given in the CM-system by the expressions 

<(2 



a.b 



x^(-p')(-i7')AV(-p) 
=:-Y,glD^\ko,k) [A'^7']®[AS^] 

a,b 

=: i[7A(-p')®^(p')]V^(p.)(fco,fc)[0(-p)®^(p)] (21) 



i-^[l(fco,fc) 



'(Pv) 
a.b 



X7/;(-p')7'7''(-»(-fc.))AbV'(-p) 



:= i[Vi(-p')®Vi(p')]^^M(fco,fe)[7A(-p)®7A(p)] (22) 

in case of pseudoscalar and pseudovector coupling, respec- 
tively. Here the meson propagator is given by 



(23) 



with fc^ — the momentum transfer. 

In instantaneous approximation we set ko = 0. From 
Eqs. |21I22] we extract the corresponding potentials in 
momentum space 



k\^+fil 



7^07^ 



(24) 



for pseudoscalar coupling and 

-1 



(25) 



= El^[A"®A'^],,„ , 

E gg [A°g)A°1 
477Z2^ ^|fc|2+^2 



(7^7 • k)(g)(j^j ■ fe) 
(7^7 • fc)(g)(7^7 • k) 



for pseudovector coupling, where k :— . 

As it stands, the expression for the potential in the 
instantaneous approximation for pseudoscalar coupling 
leads, after Fourier transformation, to a local Yukawa po- 
tential in configuration space with the usual range given 
by the mass of the exchanged pseudoscalar meson. For 
pseudovector coupling the non-relativistic approximation 
to the Fourier transform leads to the usual spin-spin con- 
tact interaction together with the usual tensor force. In the 
simplest form adopted by the Graz group [10 l^l fTH - fTT] 
the latter were ignored, in addition the contact term was 
regularised by a Gaussian function and the Yukawa terms 
were regularised to avoid singularities at the origin. 

In view of this and the instantaneous approximation 
we decided to parametrise the new flavour dependent in- 
teraction purely phenomenologically as a local potential 
in conflguration space, its simple form given by 



V<^'\x) = J2gl [A«7'®A"7'] vxA^) 



(26) 



where vx{x) is the Gaussian form given in Eq. (TH]. Other 
Dirac structures, such as [7^7 • x (8)7^7 • 2] were tried, 
but were found to be less effective. Results for the meson 
exchange form of the interaction as given by Eqs. (|24I25I) 
will be briefly discussed in sectional 
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Table 1. The hst of baryon resonances of which the masses 
were used to determine the model parameters in a least-squares 
fit where every resonance was attributed a weight reciprocal to 
its uncertainty in its position as given in 9_. Nominal masses 
are given in MeV. 



A N 


A/ E E Q 


531 (1620) Sii 


(1535) S2i(1620) Sii(1309) 


.SoinQoo^ Sii 




P3i(1750) Pii 


(939) Poi(1116) Pi3(1530) Poi(1672) 






P33(1232) 


Poi(1600) 


P33(1600) 


P2i(1660) 


P33(1920) 




7333(1700) 


Di3(1820) 


£'33(1940) 




P'35(1930) 




F35(1905) 




^35(2000) 




^37(1950) 




G37(2200) 




G39(24 00) 




-ff39(2300) 




Ha 11 (2420) 




7^315(2950) 



4 Mass spectra 
4.1 Model parameters 

The resulting baryon mass spectra were obtained by fit- 
ting the parameters of the model, viz. the offset a and 
slope h of the confinement potential, the constituent quark 
masses m„ = ruu = md and to^, the strengths of the in- 
stanton induced force, and gns as well as the strengths 
of the additional flavour dependent interaction, given by 
gs and go for flavour octet and flavour singlet exchange 
(thus assuming SU{i) symmetry) to a selection of baryon 
resonances, see table [T]. The range A given to the instan- 
ton induced force was kept to the value used in [3][S] and is 
roughly in accordance with typical instanton sizes. The op- 
timal value for the range of the additional flavour depen- 
dent interaction was found to be Ag = Aq « 0.25 fm and 
thus turned out to be of rather short range. A comparison 
of the parameters obtained with the parameters of model 
A of [11[5] is given in table [5]. The parameters need some 
comments: In the original paper [4l[5] the Dirac structures 
(i.e. the spin dependence) of the confinement potential 
were taken to be Pq = ■j(l(8)l<8il+7°(8i7°<8il + cycl. perm.) 
and Fs = ^(-l®!®! + 7°(g)7° (g) 1 -f- cycl. perm.) for the 
offset and slope, respectively, and were considered to build 
a 3-body kernel. In the present model including the ad- 
ditional octet and singlet flavour exchange potential into 
account we obtained the best results with Fq = 1®!®! 
and Pg = 7° €57'^ and treating the interaction correspond- 
ing to the latter term as a 2-body interaction. This of 
course impedes a direct comparison of the corresponding 



Table 2. Model parameters for the current model C in com- 
parison to those of model A of [H[S]. Some of the parameters 
have been slightly changed with respect to the original values 
(listed in brackets) of [US], since the calculation has been per- 
formed with higher numerical accuracy by taking more basis 
states in the diagonalisation of the Salpeter Hamiltonian, see 
also text. 



parameter 




model C 


model A 


masses 


m„ [MeV] 
rris [MeV] 


325.0 
600.0 


330.0 
670.0 


confinement 


a [MeV] 
b [MeV/fm] 


-366.78 
212.81 


-734.6 
[-744.0] 
453.6 
[440.0] 


instanton 

induced 
interaction 


ff„„ [MeV fm^] 

gns [MeV fm^] 
A [fm] 


341.49 

273.55 
0.4 


130.3 
[136.0] 
81.8 
[96.0] 
0.4 


octet 
exchange 

singlet 
exchange 


ll [MeV fm=*] 

ll [MeV fm'*] 
Ag = A() [fm] 


100.86 

1897.43 
0.25 





parameters. Furthermore it was found that the strengths 
of the instanton induced interaction is roughly tripled 
when compared to the original values. Note that the ad- 
ditional flavour exchange interaction has the same spin- 
flavour dependence as parts of the former interaction. The 
flavour singlet exchange could effectively be considered as 
an other spin dependent part of the confinement potential. 
Possibly this explains the extraordinary large coupling in 
this case. In summary it thus must be conceded that the 
present treatment is phenomenological altogether and that 
here unfortunately the relation to more fundamental QCD 
parameters, such as instanton couplings and string tension 
is lost. Nevertheless with only 10 parameters we consider 
the present treatment to be effective especially in view of 
its merits in the improved description of some resonances 
to be discussed below. 



4.2 The A and the J? spectrum 

In fig. [5] we compare the results from the present calcu- 
lation (model C) (right side of each column) with experi- 
mental data from the Particle Data Group [3] (central in 
each column) and with the results from model ^ of [3] 
(left side in each column) . The parameters used are listed 
in table [5]. The spectrum of the A (see fig.H]) and f2 (see 
right panel of fig. |6l) resonances is determined by the con- 
finement potential and the flavour exchange interaction 
only, since the instanton induced interaction does not act 
on flavour symmetric states. Concerning the positive par- 
ity resonances we see that in the present calculation we can 
now indeed account for the low position of the A3 + (1600) 
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Fig. 2. Comparison of the zi-Spectrum calculated within the present model C (right side of each column) with experimental 
data from the Particle Data Group [2] (central in each column) and with the results from model ^ of [4] (left side in each 
column), note the caption to flg. [T] Lines indicate the resonance position (mass) with the mass uncertainty represented by a 
shaded box and the rating of 9 indicated by stars. The small numbers give the mass in MeV. J and tt denote total angular 
momentum and parity, respectively. 



resonance. In addition the next excitations in this chan- 
nel now lie closer to 2000 MeV in better agreement with 
experimental data, as is also the case for the splitting of 
the two Ai+ resonances. Note, however, that these states 

2 

were included in the parameter fit. Additionally there is 
support for a parity doublet Z\3 + (1920) and Z\3-(1940) 

2 2 

as argued by pB] . 

Likewise, we can now account for the excited nega- 
tive parity resonances: Z\i- (1900) ,Z\3- (1940) and Aa- 

2 2 2 

(1930) and even find two states in the Aa- channel which 

2 

could correspond to the poorly established Z\3-(1940) 

2 

state. 

In view of the near degeneracy of the Ai- (1900), 

2 

A3- (1940) and Z\5-(1930) states it is tempting to clas- 

2 2 
sify these in a non-relativistic scheme as a total spin S =^ 

|, total quark angular momentum L = 1 multiplet, which, 
because of total isospin / = | must then belong to a 
(56, 1~) multiplet, which is lowered with respect to the 
bulk of the other negative parity states that in an oscilla- 
tor classification would be attributed to the = 3 band. 



Obviously this is not supported by the calculations: As 
table [3] shows, although the lowest J'^ = | resonance has 
a dominant component in this multiplet, the second ex- 
cited — ^ , I resonances have dominant components 
in the (70, 1~) multiplet; indeed the third excited states in 
these channels can be attributed to the (56, 1~) multiplet. 
Otherwise the description and in particular the Z\-Regge- 
trajectory are of a similar quality as in the original model 
A. 



Concerning the J7 spectrum, see fig. El (right panel), 

apart from the appearance of an excited ^23+ state at 

2 

2021 MeV no spectacular changes in the predictions with 
respect to the original model A were found. Note that the 
present model predicts that this state is almost degenerate 
with the first negative parity states fii- at 2008 MeV and 

Hs- at 1983 MeV. 

2 

In table S] we have summarised the calculation of A- 
resonances. 
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Fig. 3. Comparison of the A'^-Spectrum calculated within the present model C (right side of each column) with experimental 
data from the Particle Data Group [9] (central in each column) and with the results from model .A of [4] (left side in each 
column). See also caption to fig. [2] 



4.3 The N spectrum 

In fig. [3] we present the results for tlie nucleon spectrum. 
As was tlie case for the A spectrum, comparing to the 
results from the former model A we indeed obtain an im- 
proved description of the position of the first excited state 
with the same quantum numbers as the ground state, the 
so called Roper resonance, while at the same time improv- 
ing also on the position of the first excited negative parity 
resonances .r = (|-)^, {\~)^ and (|")^ . 

With the exception of the = | state, which compared 
to model A is shifted upwards by approximately 50 MeV, 
the description of all known excited states is of a similar 
quality as that of model A. In particular the position of 
the lowest J'^ = | resonance is still underestimated by 
more than 100 MeV. 

In the following we compare the predictions obtained 
in model C for nucleon resonances with J < | and masses 
larger than 1.8 GeV with new results obtained in the 
Bonn-Gatchina analyses as reported in [201 HI]: In par- 
ticular in [3^ a fourth J'^ = i"*" state was found, called 
A^i + (1875) which could correspond to our calculated state 

2 

at 1893 MeV. Furthermore the analysis contains two 
J"" ^ 1^ states, called N{1900)Pi3 and 7V3 + (1975), which 



might be identified with the model C states calculated at 
1825 MeV and 1945 MeV (or 1966 MeV), respectively. 
Concerning the negative parity states, in a new J"^ = 

I state was found (iVi~(1895)) which could be identi- 
fied with the calculated state at 1851 MeV (or with that at 
1881 MeV). In addition two J'^ = |" states were found: 
iV3~-(1875) could correspond to the calculated states at 

1853 MeV (or at 1934 MeV) and TVs -(2150) with one of 
the three states with calculated masses 2073 MeV, 2091 
MeV and 2137 MeV. The new TVs -(2060) state reported 

2 

in is closest to the states calculated at 1945 MeV and 
at 2007 MeV. Finally, for J'^ = |+ the analysis is am- 
biguous: Although a solution with a single pole around 
2.1 GeV is not excluded, solutions with 2 poles, either 
an ill-defined pole in the 1800-1950 MeV mass region and 
one at nearly 2.2 GeV or two close poles at approximately 
2.0 GeV were found and could correspond to the model C 
states calculated at 1930 MeV and 1983 MeV. Note that 
these new resonances were not included in the parame- 
ter fit (see table [Ij . An overview of the identification of 
nucleon resonances is given in table [5] 
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Fig. 4. Comparison of tlie yl-Spectrum calculated within the present model C (right side of each column) with experimental 
data from the Particle Data Group [9] (central in each column) and with the results from model ^ of [4] (left side in each 
column), see also caption to fig. [D 



4.4 Hyperon spectra 

The resulting spectra for hyperon resonances, viz. the A ,E 
and E states are depicted in figs. |H [S] and 0, respec- 
tively. Again we indeed find an improved description of the 

"Roper-like" resonances (1600) and S3 + (1660). Note, 

2 2 

however, that both were used to determine the model pa- 
rameters. Concerning the negative parity resonances, al- 
though we do find an acceptable description of the A reso- 
nances with J'^ = I , I and | , also the new calculat- 
ion can not account for the low position of the Ai^ (1405) 

2 

resonance, which now is 200 MeV below the calculated po- 
sition. In our opinion this underlines the conclusion, that 
this state cannot indeed be accounted for in terms of a 
excitation alone and that its position is determined by a 
strong coupling of a "bare" q'^ state to meson-baryon de- 
cay channels due to the proximity of the ^A^-threshold, 
see also ref. [22H24] for a description of this state in a 
chiral unitary approach. 

Concerning the S resonances, with respect to model 
of [S] mainly the prediction for the excited state 
at 1765 MeV is 100 MeV lower. To a lesser extend this 
also holds for the excited ^(1 ) which is now predicted 
at 1889 MeV. 



5 Electromagnetic properties 

As has been elaborated in 6.^ in lowest order the transition 
current matrix element for an initial baryon state with 
four-momentum P' — M = [M, 0) in its rest frame and a 
final baryon state with four-momentum P is given by the 
expression 

(F|,-(0)|M) =-^J§^J§^.mPe,Pn~h) 

r^APi,Pv)- (27) 

Here q denotes the momentum transfer, q is the charge 
operator and 

<Pij{p'^,p'^) (28) 

is the vertex function in the rest frame of the baryon 
with Va the projection of the instantaneous interaction 
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Fig. 5. Comparison of the X'-Spectrum calculated within the present model C (right side of each column) with experimental 
data from the Particle Data Group [9] (central in each column) and with the results from model ^ of [4] (left side in each 
column), see also caption to fig. [D 



kernels onto the subspace of purely positive and purely 
negative energy components only, see in particular Ap- 
pendix A of jB^ for details. The vertex for a general four- 
momentum on the mass shell can be obtained by an ap- 
propriate Lorentz-boost. 

Accordingly, the Sachs form factors are given by 



^ {N,P,^\j^mN,M,^) 
'^mW ) — 



(29a) 
(29b) 



where | A'', P, A) denotes a nucleon state with four-momen- 
tum P and helicity A. Furthermore — —q^ and := 



jf ±ijf as well as |P|2 = Q2 
the axial vector form factor is given by 



4M2 



. Likewise, 



{p,P,^\4{0)\n,M,~^) 



(30) 



where again :— jj^± ij^ with j'^ the axial current 
operator, whose matrix elements are given by Eq. (j27p 
after the formal substitution q . Of course, the 



normalisations of the form factors is such that the static 
magnetic moments and the axial coupling are given by 



Gm(O' = 0), ffA :=G^(0' =0). 



5.1 The electric form factors of the nucleon 



In fig. [7] and [5] we display the electric proton and neutron 
form factor, respectively, up to a momentum transfer of 
g2 = 6.0 GeV^ The black solid curve is the result of the 
present model C, the blue dashed curve is the result ob- 
tained with the parameters of model as in [5], albeit 
with a better numerical precision, see the end of Subsec- 
tion [ 



Although the electric form factor of the proton, see 
fig. [71 as calculated with model .4 in [B] fell too steeply 
in comparison to experimental data, with the present in- 
teraction wc find a much improved shape which yields a 
satisfactory description even up to momentum transfers 
of 6 GcV^ . Indeed, in contrast to model A, which mainly 
failed with respect to the isovector part of the form factor, 
in the present model C this form factor shows an almost 
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Fig. 6. Comparison of the H-Spectrum (flrst eight columns) and the J7-Spectrum (rightmost eight columns) calculated within 
the present model C (right side of each column) with the experimental data from the Particle Data Group |9j (central in each 
column) and with the results from model A of ,4, (left side in each column), see also caption to fig. [D 



perfect dipole shape with the parametrisation 

(1 + Q2/M2)2' (31) 

taken from [25l[26] with = 0.71 GeV^ 

The resulting electric neutron form factor, see fig. HI 
has a maximum at approximately the experimental value 
of « 0.4 GeV^ but underestimates the experimental 
data from [SI] by about the same amount as the ear- 
lier calculation overestimated the data. However, the pre- 
diction of model C is very similar to the predictions of 
the Graz group _17J and [72 for the Goldstone-boson- 
exchange quark models. The corresponding charge radii 
are given in table [H]. As for the form factor the result- 
ing squared charge radius of the neutron is calculated too 
small by a factor of two. Also the r.m.s. proton radius is 
slightly smaller than the experimental value. 

In fig. inland [TU]we display the magnetic proton- and 
neutron form factor up to a momentum transfer of — 
6.0 GeV^ , respectively. Again, the black solid curve is the 
result of the present model C, the blue dashed curve is the 
result obtained with the parameters of model A, as in [5] , 
albeit with a better numerical precision, see the end of 
Subsection 12.21 . 

Whereas in the original calculation (model A of [Bj) 
the absolute value of these form factors dropped slightly 



too fast as a function of the momentum transfer, in the 
present calculation we now find a very good description 
even at the highest momentum transfers. Only at low mo- 
mentum transfer the values are too small as is reflected by 
the rather small values for the various magnetic radii, see 
table[B]and the too small values of the calculated magnetic 
moments. Note, however that the ratio fip/iin ~ 1.597 
for model C slightly changes (previously fip/nn ~ 1.605 
for model A) and is slightly larger than the experimental 
value fip/ fin ~ 1.46 ; all values are remarkably close to the 
non-relativistic constituent quark model value fip/fin = 
I . The magnetic moments of fiavour octet and decuplet 
baryons has been calculated accordingly to the method 
outlined in [8,. The results are compared to experimental 
values in Table [7] and |S1 respectively. As a consequence of 
the better description of the momentum transfer depen- 
dencies in the individual form factors, we now also find 
an improved description of the momentum transfer de- 
pendence of the form factor ratio fip G^^/G\j{Q'^) , which 
has been the focus on the discussion whether two-photon 
amplitudes are relevant for the discrepancy [55] found be- 
tween recent measurements based on polarisation data 
(red data points of fig. [TT|) [50 5^ [51HBU] versus the tradi- 
tional Rosenbluth separation (black data points of fig. fTTj) . 
see e.g. [4ni45M7] . Whereas in the original model A this 
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Table 3. Muitiplet decomposition of amplitudes of negative 
parity Zi-resonances. For eacli amplitude the contribution to 
the Salpeter norm, see is given in % , in each row the upper 
line and the lower line give the positive and negative energy 
contribution, respectively. States are labeled by the calculated 
mass and denotes total angular momentum and parity, 
label amplitudes with spin S, flavour representa- 
tion with dimension SU{G) representation with dimension 
D . The dominant contribution is underlined. 



Table 5. Comparison of experimental [9] and calculated 
masses in MeV of A''-resonances. The corresponding spectra 
are shown in fig. |3]. 



J" 


Mass 


pos. 


Ho J [56] 


Ho./ [70] 






neg. 


Ho./ [56] 


Ho J [70] 


1 - 

2 


1635 


98.9 


6.6 


92.3 






1.1 


0.8 


0.4 


1 - 
2 


1932 


98.5 


14.2 


84.3 






1.5 


1.0 


0.5 


1 - 
2 


2041 


99.0 


81.6 


17.3 






1.0 


0.5 


0.5 


3 - 
2 


1592 


98.2 


9.9 


88.4 






1.8 


0.9 


0.9 


3 - 
2 


1871 


97.9 


24.5 


73.3 






2.1 


1.4 


0.7 


3 - 
2 


1944 


98.6 


65.1 


33.5 






1.4 


0.9 


0.6 


5 - 
2 


2013 


98.8 


89.3 


9.6 






1.2 


0.6 


0.6 


5 - 
2 


2136 


99.2 


16.8 


82.4 






0.8 


0.3 


0.5 


5 - 
2 


2166 


99.1 


86.6 


12.4 






1.0 


0.4 


0.6 



Table 4. Comparison of experimental i9j and calculated 
masses in MeV of zl-resonances. The corresponding spectra 
are shown in fig. [2]. 



exp. 


rating model C 


exp. rating 


model C 


5'3i(1620) 
S'3i(2150) 


1634 
* 2040/2109 


531(1900) 


1932 


P3i(1750) 
^33(1232) 
P33(1920) 


1653 
1231 

*" 1894/1925 


P3i(1910) 

P33(1600) 


1888 
1559 


D33(1700) 
D35(1930) 


1591 
2012 


L»33(1940) 


1871/1944 


F35(1905) 
^37(1950) 


1897 
1936 


^35(20 00) 
f 37 (2390) 


1961 
many res. 


G37(2200) 


* 2126/2198 


G39(2400) 


2232 


^39(23 00) 


2319^3,11(2420) 


2372 


73,13(2750) 


2584 






/S'3,i5(2950) 


" 2724/2909 







ratio fell much too steep, we now find in model C a much, 
better description of this quantity, see fig. [TT] for a com- 



exp. 


rating model C exp. rating 


model C 


Sii(1535) 
5'ii(2090) 


1484 Sii(1650) 
* many res. 


1672 


Pii(939) 
Ai(1710) 
As (1720) 


945 Pii(1440) 
1709 Pii(2100) 
1703 Pi3(1900) 


1440 
many res. 
1825 


Di3(1520) 
Di3(2080) 
Di5(2200) 


1534 Dig (1700) 
many res. 1)15(1675) 
** many res. 


1685 
1667 


Fi5(1680) 
Fi7(1990) 


1761 Pi5(2000) 
2001 


1930/1983 


Gi7(2190) 


2011 Gi9(2250) 


2165 


ffi9(2220) 


2192 




/i, 11 (2600) 


2377 




7^1,13(2700) 


2543 






3 4 

[GeV^] 

Fig. 7. The electric form factor of the proton divided by 
the dipole form Gd(Q^), Eq. MMD-Data are taken from 

Mergell et al. supplemented by data from Christy et 

al. [27] and Qattan et al. [28] . The solid black line represents 
the results from the present model C; the dashed blue line those 
from model ^ of [6], albeit recalculated with higher numerical 
precision. Red data points are taken from polarisation experi- 
ments and black ones are obtained by Rosenbluth separation. 



parison with various data. Up to « 3 GeV^ we indeed 
find the observed linear dependence. 

Finally, the axial form factor, see fig. [T^l was already 
very well described in model ^ of 6 . Although falling 
slightly less steeply, the present calculation still gives a 
very satisfactory description of the data also at higher 
momentum transfers in the same manner as in [73H75j . As 
for the magnetic moments the value of the axial coupling 
constant is too small, but of course much better than the 
non-relativistic constituent quark model result g^i = | . 
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Fig. 8. Tlie electric form factor of the neutron. MMD-Data 
are taken from the compilation of Mergell et al. [25]. The solid 
black line represents the results from the present model C; the 
green- and brown-dashed lines corresponds to the Yukawa-like 
models with pseudoscalar and/or pseudovector coupling, see 
Eqs. (|24|l and (|25|l respectively; the dashed blue line in the 
result from model ^ of [B], albeit recalculated with higher nu- 
merical precision. Red data points are taken from polarisation 
experiments and black ones are obtained by Rosenbluth sepa- 
ration. 

' ' ' Wimp [25] i i i 
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0. Morten A i£ 

1.2 - 1 • model C " 
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[GeV^] 

Fig. 9. The magnetic form factor of the proton divided by 
the dipole form Gd{Q^), Eq. pi|l and the magnetic moment 
of the proton jj,p = 2.793 /xjv. MMD-Data are taken from the 
compilation of Mergell et al. Additionally, polarisation 

experiments are marked by red. The black marked data points 
are obtained by Rosenbluth separation. 

The axial form factor, presented in fig. [T2j is divided by 
the axial dipole form 

'^^w'' = (i + oW - 

with the parameters = 1.014 ± 0.014 GeV and = 
1.267 taken from Bodek et al. [26,. 

In summary we find, that the new model C, apart from 
some improvements in the description of the excitation 



Table 6. Static properties of the nucleon. The values in paren- 
theses are as reported in 0, the values on top of these are 
obtained within the same model A but with higher numerical 
accuracy. The columns PS and PV are the results obtained 
with the additional interaction kernels according to Eqs. (|24p 
and (|25|) . respectively, see also Section [6] for a brief remark. 
The static values are extrapolated from a dipole shape-like fit. 





Model A 


PS 


PV 


model C 


exp. 


ref. 


Mp[Miv] 


2.76 
[2.74] 


2.49 


2.39 


2.54 


2.793 


[5] 




-1.71 
[-1.70] 


-1.59 


-1.54 


-1.59 


-1.913 


lej 


V(r2)|[fm] 


0.91 
[0.82] 


0.83 


0.68 


0.81 


0.847 


.25^ 




-0.20 
[-0.11] 


0.01 


0.08 


-0.06 


-0.123 
±0.004 






0.90 
[0.91] 


0.81 


0.68 


0.78 


0.836 


m 




0.84 
[0.86] 


0.79 


0.67 


0.75 


0.889 


m 


9A 


1.22 
[1.21] 


1.17 


1.14 


1.13 


1.267 
±0.0035 




V(r2)A[fm] 


0.68 
[0.62] 


0.64 


0.48 


0.57 


0.67 
±0.01 





Table 7. Octet hyperon magnetic moments /i for model A and 
C calculated as in [S] . The values are given in units of ^jv . 



hyperon 


model A 


model C 


PDG 9 


A 


-0.606 


-0.577 


-0.613 ±0.004 


E+ 


2.510 


2.309 


2.458 ±0.010 




0.743 


0.701 




s- 


-1.013 


-0.908 


-1.160 ±0.025 


=■0 


-1.324 


-1.240 


-1.250 ±0.014 




-0.533 


-0.532 


-0.651±0.0025 



spectra at the expense of additional parameters of a phe- 
nomenologically introduced flavour dependent interaction 
does allow for a parameter-free description of electromag- 
netic ground state properties of a similar overall quality as 
has been obtained before, with some distinctive improve- 
ments on the momentum transfer dependence of various 
form factors. A discussion of the momentum dependence 
of helicity amplitudes for the electromagnetic excitation 
of baryon resonances will be given in a subsequent pa- 
per [75] . 

6 Summary and conclusion 

In the present paper we have tried to demonstrate that by 
introducing an additional flavour dependent interaction. 
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Table 8. Decuplet hyperon magnetic moments /i for model A 
and C calculated as in [S]. The values are given in units of /ijv. 



hyperon 


model A 


model C 


PDG ^9, 




4.241 


4.238 


3.7 to 7.5 


A -t- 


2.121 


2.119 


2.7ti;3 ± 1-5 ± 3 


A 


0.0 


0.0 




A- 


-2.121 


-2.119 




E*+ 


2.567 


2.431 






0.275 


0.205 




s*- 


-2.017 


-2.021 






0.607 


0.474 






-1.865 


-1.765 




n- 


-1.675 


-1.577 


-2.02 ± 0.05 




Fig. 10. The magnetic form factor of the neutron divided by 
the dipole form Gd{Q^), Eq. (|3ip and the magnetic moment of 
the neutron /i„ = —1.913 pijv. MMD-Data are taken from the 
compilation by Mergell et al. |25) and from more recent results 
from MAMI [42144) . Additionally, polarisation experiments are 
marked by red data points. The black marked ones are obtained 
by Rosenbluth separation. 



parametrised with a Gaussian radial dependence with an 
universal range and two couplings for flavour octet and 
flavour singlet exchange, it is possible to improve upon 
some deficiencies found in a former relativistically covari- 
ant constituent quark model treatment of baryonic exci- 
tation spectra based on (an instantaneous formulation of) 
the Bethe-Salpeter equation. These improvements include: 

— A better description of excited negative parity states 
slightly below 2 GeV in the A spectrum; 

— A better description of the position of the first scalar, 
isoscalar excitation of the ground state in all light-fla- 
vour sectors; 

— An improved description of the momentum depend- 
ence of electromagnetic form factors of ground states 
without the introduction of any additional parameters. 



s 



3, 




Fig. 11. The ratio /Xp G^/G^ compared to recent JLAB data 
(see legend). In the insert the low momentum transfer region is 
enlarged. The solid black line is the present result, the dashed 
blue line the result in model A . Red data points are taken from 
polarisation experiments and the black ones are obtained from 
Rosenbluth separation. 




Amaldi 
Brauel 
Bloom 
Del Guerra 
Joos 
Baker 
Miller 
Kitagaki 83 
Kitagaki 90 
Allasia 
model ^ 



61 
62 
63 

13 

65 
66 
67 
68 
69 
TO 
6 



model C 



[GeV^ 



Fig. 12. The axial form factor of the nucleon divided by the 
axial dipole form in Eq. p2[) and the axial coupling gA = 
— 1.267. The black solid line is the present result, the blue 
dashed line the result in model A . Experimental data are taken 
from the compilation by Bernard et al. [7T] . 



It must be conceded that this additional interaction was 
introduced purely phenomenologically and required a 
drastic modification of the parametrisation of confinement 
and the other flavour dependent interaction of the origi- 
nal model, which had a form as inferred from instanton 
effects. In spite of this, with only 10 parameters in total 
we still consider this to be an effective description of the 
multitude of resonances found for baryons made out of 
light flavoured quarks. 

Nevertheless, it would have been preferred, if the ad- 
ditional flavour dependent interaction could be related to 
a genuine physical process, such as light pseudoscalar me- 
son exchange. Indeed, we tried to find parametrisations 
of confinement and parameters of the instanton induced 
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Fig. 13. Comparison of the Zi-Spectrum (first six columns) and the A''-Spectrum (rightmost six columns) calculated within 
the PS and PV coupled models from Eqs. H24[l and H25|) (right and left side of each column) and with the experimental data 
from the Particle Data Group (central in each column). 



interaction, that could be combined with interaction ker- 
nels as given by the expressions in Eq. ([M)) or How- 
ever, for pseudoscalar coupling (PS) of a meson nonet to 
the quarks, we could find a description of the mass spec- 
tra with similar features and of a similar quality as in the 
present model C only at the expense of introducing flavour 
SU(3) symmetry breaking and thus introducing more pa- 
rameters. The latter was also found to be the case for 
pseudovector coupling (PV) where, moreover, no signifi- 
cant improvement concerning the deficiencies in the spec- 
tra mentioned above was found. For the sake of complete- 
ness, the nucleon spectrum for PS and PV coupled models 
from Eqs. (IMl) and is shown in fig. [T31 Furthermore, 
with both Ansatze we were not able to reproduce in partic- 
ular the electric neutron form factor, see fig. [5] and table [5] 
for some typical results. Accordingly, we dismissed these 
possibilities and preferred the phenomenological approach 
of model C discussed in the paper. 



A parameter- free calculation of longitudinal and trans- 
verse helicity amplitudes for electro-excitation is presently 
performed and the results will be discussed in a subsequent 
paper [7B] . 



Acknowledgments 

Stimulating discussions with E. Klempt, A. V. Sarant- 
sev and U. Thoma within the framework of the DFG 
supported Collaborative Research Centre SFB/TR16 are 
gratefully acknowledged. 



References 

1. R.G. Edwards, J.J. Dudek, D.G. Richards, S.J. Wallace, 
JLAB-THY-1370 arXiv:1104._51 52vl] (2011) 

2. Huey-Wen Lin, NT@UW-ll-09 |arXiv:1106.T608V l] 
(2011) 

3. U. Loring, K. Kretzschmar, B.C. Metsch, H.R. Petry, Eur. 
Phys. J. A 10, 309 (2001). 

4. U. Loring, B.C. Metsch, H.R. Petry, Eur. Phys. J. A 10, 
395 (2001). 

5. U. Loring, B.C. Metsch, H.R. Petry, Eur. Phys. J. A 10, 
447 (2001). 

6. D. Merten, U. Loring, K. Kretzschmar, B.C. Metsch, H.R. 
Petry, Eur. Phys. J. A 14, 477 (2002). 

7. T. van Cauteren, D. Merten, T. Corthals, S. Janssen, B.C. 
Metsch and H.R. Petry, Eur. Phys. J. A 20, 283 (2004). 

8. C. Haupt, B.C. Metsch, H.R. Petry, Eur. Phys. J. A 28, 
213 (2006). 



M. Ronniger et al.: Effects of a spin-flavour dependent interaction on the baryon mass spectrum 



15 



9. K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 
075021 (2010). 

10. L.Ya. Glozman, D.O. Riska, Phys. Rep. 268, 263 (1996). 

11. L.Ya. Glozman, Z. Papp, W. Plessas, K. Varga and R.F. 
Wagenbrunn, Nucl. Phys. A 623, 90 (1997). 

12. L.Ya. Glozman, Z. Papp, W. Plessas, K. Varga and R.F. 
Wagengrunn, Phys. Rev. C 57, 3406 (1998). 

13. L.Ya. Glozman, W. Plessas, K. Varga and R.F. Wagen- 
grunn, Phys. Rev. D 58, 094030 (1998). 

14. L. TheuBl, R.F. Wagenbrunn, B. Desplanques and W. 
Plessas, Eur. Phys. J. A 12, 91 (2001). 

15. K. Glantschnig, R. Kainhofer, W. Plessas , B. Sengl and 
R. F. Wagenbrunn, Eur. Phys. J. A 23, 507 (2005). 

16. T. Melde, W. Plessas, and B. Sengl, Phys. Rev. D 77, 
114002 (2008). 

17. W. Plessas and T. Melde, AIP Conf. Proc. 1056, 15 
(2008). 

18. G. Caia, J.W. Durso, Ch. Elster, J. Haidenbauer, A. Sibirt- 
sev and J. Speth, Phys. Rev. C 66, 044006 (2002). 

19. L Horn et al, Phys. Rev. Lett. 101, 202002 (2008). 

20. A.V. Anisovich, E. Klempt, V.A. Nikonov, A.V. Sarantsev 
and U. Thoma, Eur. Phys. J. A 47, 27 (2011). 

21. A.V. Anisovich, V.A. Nikonov, A.V. Sarantsev, U. Thoma 
and E. Klempt, arXiv:1109.0970vl] (2011) 

22. D. Jido, J.A. Oiler, E. Oset, A. Ramos, U. -G. Meifiner, 
Nucl. Phys. A 725, 181 (2003). 

23. T. Hyodo, D. Jido and L. Roca, Phys. Rev. D 77, 056010 
(2008). 

24. T. Hyodo, D. Jido and A. Hosaka, Phys. Rev. C 78, 025203 
(2008). 

25. P. Mergell, U.-G. Meifiner and D. Drechsel, Nucl. Phys. A 
596, 367 (1996). 

26. A. Bodek, S. Avvakumov, R. Bradford, and H. Budd, J. 
Phys. Conf. Ser. 110, 082004 (2008). 

27. M.E. Christy et al, Phys. Rev. C 70, 015206 (2004). 

28. LA. Qattan et al, Phys. Rev. Lett. 94, 142301 (2005). 

29. T. Eden et al, Phys. Rev. C 50, 1749 (1994). 

30. C. Herberg et al, Eur. Phys. J. A 5, 131 (1999). 

31. M. Ostrick et al, Phys. Rev. Lett. 83, 276 (1999). 

32. L Passchier et al, Phys. Rev. Lett. 82, 4988 (1999). 

33. D. Rohe et al, Phys. Rev. Lett. 83, 21 (1999). 

34. R. Schiavilla et al, Phys. Rev. C 64, 041002 (2001). 

35. J. Golak et al, Phys. Rev. C 63, 034006 (2001). 

36. H. Zhu et al, Phys. Rev. Lett. 87, 081801 (2001). 

37. R. Madey et al, Phys. Rev. Lett. 91, 122002 (2003). 

38. G. Warren et al, Phys. Rev. Lett. 92, 042301 (2004). 

39. D.L Glazier et al, Eur. Phys. J. A 24, 101 (2005). 

40. R. Alarcon et al, Eur. Phys. J. A 31, 588 (2007). 

41. T. Bartel et al, Nucl. Phys. B 58, 469 (1973). 

42. H. Anklin et al, Phys. Lett. B 428, 248 (1998). 

43. W. Xu et al, Phys. Rev. Lett. 85 (2000). 

44. G. Kubon et al, Phys. Lett. B 524, 26 (2002). 

45. L.E. Price et al, Phys. Rev. D 4 (1971). 

46. Ch. Berger , V. Burkert , G. Knop , B. Langenbeck and 
K. Rith, Phys. Lett. B 35, 1 (1971). 

47. R.C. Walker et al, Phys. Rev. D 49 (1994). 

48. L. Andivahis et al, Phys. Rev. D 50 (1994). 

49. M. Vanderhaeghen, Nucl. Phys. A 755, 269c (2005). 

50. B.D. Milbrath et al, Phys. Rev. Lett. 80 (1998). 

51. M.K. Jones et al, Phys. Rev. Lett. 84 (2000). 

52. O. Gayou et al, Phys. Rev. C 64, 038202 (2001). 

53. T. Pospischil et al, Eur. Phys. J. A 12, 125 (2001). 

54. O. Gayou et al, Phys. Rev. Lett. 88 (2002). 



55. V. Punjabi et al, Phys. Rev. C 71, 055202 (2005). 

56. B. Hu et al, Phys. Rev. C 73, 064004 (2006). 

57. C.B. Crawford et al, Phys. Rev. Lett. 98, 052301 (2007). 

58. D.W. Higinbotham, AIP Conf. Proc. 1257, 637 (2010). 

59. G. Ron et al, arXiv:1103.5784vl [nucl-ex] {2QII). 

60. X. Zhan et al, JLAB-PHY-11-1311, 5 (2011). 

61. E. Amaldi et al, Phys. Lett. B 41, 216 (1972). 

62. P. Brauel et al, Phys. Lett. B 45, 389 (1973). 

63. E.D. Bloom et al, Phys. Rev. Lett. 30, 1186 (1973). 

64. A. Del Guerra et al, Nucl. Phys. B 99, 253 (1975). 

65. P. Joos et al, Phys. Lett. B 62, 230 (1976). 

66. N.J. Baker et al, Phys. Rev. D 23, 2499 (1981). 

67. K.L. Miller et al, Phys. Rev. D 26, 537 (1982). 

68. T. Kitagaki et al, Phys. Rev. D 28, 436 (1983). 

69. T. Kitagaki et al, Phys. Rev. D 42, 1331 (1990). 

70. D. AUasia et al, Nucl. Phys. B 343, 285 (1990). 

71. V. Bernard, L. Elouadrhiri, and U.-G. Meifiner, J. Phys. 
G 28, 1 (2002). 

72. T. Melde, K. Berger, L. Canton, W. Plessas and R.F. Wa- 
genbrunn, Phys. Rev. D 76, 074020 (2007). 

73. L.Ya. Glozman, M. Radici, R.F. Wagenbrunn, S. Boffi, W. 
Klink, W. Plessas, Phys. Lett. B 516, 183 (2001). 

74. R.F. Wagenbrunn, S. Bofli, L.Ya. Glozman, W. Klink, W. 
Plessas and M. Radici, AIP Conf. Proc. 603, 319 (2001). 

75. R.F. Wagenbrunn, S. Bofli, L.Ya. Glozman, W. Klink, W. 
Plessas, and M. Radici, Eur. Phys. J. A 18, 155 (2003). 

76. M. Ronniger, B. Ch. Metsch, in preparation. 



